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ABSTRACT 

We study the scenario of thermal leptogenesis in which the leptonic asymmetries are resonantly 
enhanced through the mixing of nearly degenerate heavy Majorana neutrinos that have mass 
differences comparable to their decay widths. Field-theoretic issues arising from the proper sub- 
traction of real intermediate states from the lepton-number-violating scattering processes are 
addressed in connection with an earlier developed resummation approach to unstable particle 
mixing in decay amplitudes. The pertinent Boltzmann equations are numerically solved after 
the enhanced heavy-neutrino self-energy effects on scatterings and the dominant gauge-mediated 
collision terms are included. We show that resonant leptogenesis can be realized with heavy 
Majorana neutrinos even as light as ~ 1 TeV, in complete accordance with the current solar and 
atmospheric neutrino data. 

PACS numbers: ll.30.Er, 14.60.St, 98.80.Cq 
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1 Introduction 



The recent results from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite 
have dramatically improved the accuracy of many cosmological parameters [1], thus sig- 
nalling a new era of precision cosmology. For the first time, the baryon-to-photon ratio of 
number densities rjB has been measured to the unprecedented precision of less than 10%. 
The reported value for rjB is [1] 

VB ^ 6.1 X 10-^° , (1.1) 

where ub — rib — ni and are the number densities of the net baryon number B and 
photons at the present epoch, respectively. 

Many theoretical models have been suggested in the literature [2,3] in order to explain 
the presently small but non-zero value of ry^ that quantifies the so-called cosmological 
Baryon Asymmetry in the Universe (BAU). One of the most attractive as well as field- 
theoretically consistent scenarios of baryogenesis is the one proposed by Fukugita and 
Yanagida [4]. In this model, out-of-equilibrium L- violating decays of singlet neutrinos Ni 
with Majorana masses considerably larger than the critical temperature Tf. 100-200 GcV 
produce initially an excess in the lepton number L. This excess in L is then converted 
into the observed B asymmetry through {B + L)-violating sphalcron interactions [5,6], 
which are in thermal equilibrium for temperatures ranging from Tc up to 10^^ GeV [7-9]. 
Many studies have been devoted to analyze in detail this scenario of baryogenesis through 
leptogenesis [10-33]. 

In the last few years, the on-going neutrino experiments, mainly at Super-K [34] and 
SNO [35], have been able to address another important question in particle and astro- 
particle physics [36]. Their analyses have offered overwhelming support to the theoretical 
idea that the ordinary neutrinos have tiny but non-zero masses and mixings [37], thereby en- 
abling them to oscillate from one type of lepton to another [38]. In the Standard Model (SM) 
neutrinos arc strictly massless. An economical as well as natural solution to this problem 
can be achieved by augmenting the SM field content with right-handed (singlet) neutrinos. 
By the same token, bare Majorana masses that violate the lepton number by two units 
arc allowed to be added to the Lagrangian. The scale of these singlet masses is rather 
model-dependent and may range from about 1 TeV in Left-Right Symmetric [39,40] or 
certain Eg [41] models up to 10^^ GeV in typical Grand Unified Theories (GUTs) such 
as SO(IO) [42,43] models. From the low-energy point of view, the large Majorana masses 
present in the complete neutrino mass matrix give rise to a kind of a seesaw mechanism [44] , 
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Figure 1: Feynman diagrams contributing to the L-violating decays of heavy Majorana 
neutrinos, Ni — > L'-'^^ , where L and $ represent lepton and Higgs-boson iso-doublets, 
respectively: (a) tree-level graph, and one-loop (b) self-energy and (c) vertex graphs. 



through which the phenomenologically favoured values for neutrino masses of order 0.1 eV 
and smaller can be explained without unnaturally suppressing the Yukawa couplings of the 
theory. 

One of the central questions that several articles have been addressing recently is to 
which extent the afore-mentioned heavy Majorana neutrinos can be responsible for both 
the observed BAU and the neutrino oscillation data, including possible data from other 
non-accelerator experiments. In this context, it has been found [19,20] that if the heavy 
singlet neutrinos have an hierarchical mass spectrum, a lower bound of about 10^-10^ GeV 
on the leptogenesis scale can be derived. In the derivation of this lower bound, the size of 
the leptonic asymmetry between the heavy Majorana neutrino decay into a lepton doublet 
L and a Higgs doublet Ni — > L$, and its respective charge and parity (CP) conjugate 
mode, Ni —>■ L^^^ , plays a key role. In other words, the larger the leptonic CP asymmetry, 
the smaller the lower bound on the leptogenesis scale becomes. 

As is shown in Fig. ^ in a Feynman-diagrammatic way, there are two one- loop graphs 
that contribute to the CP-violating leptonic asymmetry. In particular, the interference of 
the tree-level decay amplitude with the absorptive parts of the one-loop self-energy and 
vertex graphs violates CP and hence gives rise to a non-vanishing leptonic asymmetry. 
These self-energy and vertex contributions are often termed in the literature [32, 45] e- 
and e'-types of CP violation, respectively. Unlike e'-type [4,10,11], e-type CP violation 
can be considerably enhanced [12-14] through the mixing of two nearly degenerate heavy 
Majorana neutrinos. 

The fact that e-type CP violation can become several orders of magnitude larger than 
e'-tjpe CP violation might raise concerns on the validity of perturbation theory. Indeed, 
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finite-order perturbation theory breaks down if two lieavy Majorana neutrinos become 
degenerate. However, based on a field-tlieoretic approacli tliat consistently resums all the 
higher-order self-energy-enhanced diagrams, it has been shown in [14] that the Icptonic CP 
asymmetry is not only analytically well-behaved, but it can also be of order unity if two 
of the heavy Majorana neutrinos have mass differences comparable to their decay widths. 
Because of this resonant enhancement of the leptonic asymmetries, we call this scenario of 
leptogenesis resonant leptogenesis. 

An immediate consequence of resonant leptogenesis is that the singlet mass scale 
can be drastically lowered to TeV energies [14,32]. However, these previous studies have 
not considered possible limits that may arise from the presently better constrained light- 
neutrino sector. In this paper we will analyze the scenario of resonant leptogenesis in light 
of the current solar and atmospheric neutrino data [46,47]. In particular, we will show that 
resonant leptogenesis can occur with heavy Majorana neutrinos even as light as ~ 1 TeV, 
within the framework of light-neutrino scenarios with normal or inverted mass hierarchy 
and large I'e-^fj. and i'fj,-i'T mixings, namely within schemes currently suggested by neutrino 
oscillation data. 

Our predictions for the BAU are obtained after numerically solving a network of 
Boltzmann Equations (BEs) related to leptogenesis. In our analysis, we include the dom- 
inant collision terms that account for 2^2 scatterings involving the SU(2)l and U(l)y 
gauge bosons. Furthermore, the resonantly enhanced CP- violating as well as CP-conserving 
effects on the scattering processes thanks to heavy-neutrino mixing are taken into account. 
To the best of our knowledge, these two important contributions to the BEs have not been 
considered in the existing literature before. 

The proper description of the dynamics of unstable particles and their mixing phe- 
nomena is a subtle issue within the context of a field theory. To deal with this problem, 
one is compelled to rely on resummation approaches to unstable particles that consistently 
maintain all desirable field-theoretic properties, such as gauge-invariance, analyticity and 
unitarity [48,49]. In this context, a resummation approach to unstable particle mixing 
in decay amplitudes was developed in [14] which preserves CPT invariance and unitar- 
ity [32,50]. 

In this paper we address another important issue related to the proper subtraction of 
the so-called real intermediate states (RIS) from the L- violating 2-^2 scattering processes 
that result from the exchange of unstable particles in the s-channel. Such a subtraction is 
necessary in order to avoid double-counting in the BE's [51] from the already considered 
1 — > 2 decays and 2 — > 1 inverse decays of the unstable particles, namely those associated 
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with heavy Majorana neutrino decays. By examining the analytic properties of the pole and 
residue structures [49,50] of a resonant L-violating scattering amplitude, we can identify 
the part of the 2 — > 2 amplitude that contains RIS contributions only. We find that the so- 
derived resonant amplitude exhibits the very same analytic form with the one obtained with 
an earlier proposed resummation method [14]. Since the present derivation does not rely 
on resorting to a kind of Lehmann-Symanzik-Zimmermann (LSZ) reduction formalism [52] 
for the decaying unstable particle, it offers therefore a firm and independent support to the 
earlier treatment presented in [14]. 

The paper is organized as follows: in Section El we discuss the generic structure of a 
heavy Majorana- neutrino model that possesses a low singlet scale and predicts nearly de- 
generate heavy Majorana neutrinos. Employing the Froggatt-Nielsen (FN) mechanism [53], 
we also put forward a generic texture for the light neutrino mass matrices that enable an 
adequate description of the present solar and atmospheric neutrino data. In Section El 
we address field-theoretic issues that arise from the proper subtraction of RIS from the 
L-violating scattering processes. In particular, we explicitly demonstrate how the resonant 
part of the scattering amplitude is intimately related to the resummed decay amplitude 
derived earlier by means of an LSZ-type resummation approach [14]. Analytic formulae re- 
lated to the general case of three heavy-Majorana-neutrino mixing are given in Appendix 1X1 
In Section |31 we derive the relevant network of BE's for resonant leptogenesis, where the 
gauge-mediated collision terms and the resonantly enhanced CP-violating as well as CP- 
conserving contributions to scatterings due to heavy neutrino mixing are taken into account. 
Analytic expressions of reduced cross-sections for all relevant 2 — 2 scattering reactions 
are presented in Appendix El Our conclusions are summarized in Section O 



2 Low-Scale Heavy Majorana-Neutrino Model 
and Neutrino Data 

In this section, we first set up our conventions by briefly reviewing the low-energy struc- 
ture of a minimally extended SM that includes heavy Majorana neutrinos. We then put 
forward a generic scenario that predicts nearly degenerate heavy Majorana neutrinos at 
the TeV scale and can naturally be realized by means of the FN mechanism [53]. In 
this generic scenario, the light-neutrino sector admits the Large Mixing Angle (LMA) 
Mikheyev-Smirnov-Wolfenstein (MSW) [38] solution and so may explain the solar neu- 
trino data through a large Ug-i^fi mixing. The light-neutrino sector also allows for a large 
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z/^-z/r-iiiixing to account for the atmospheric neutrino anomaly. Another property of our 
generic scenario is that it leads to a mass spectrum for the light neutrinos, denoted as 
^^1,2,3 (with the mass convention rriiy^ < rriy^ < rriij,^), with normal or inverted hierarchy, 
depending on whether the lightest physical neutrino z/i has predominantly a z/g or a Ut 
component. In particular, the generic scenario can accommodate the phenomenologically 
favoured neutrino-mass differences [46,47]: 

1.4x10-^ < Am^t^ [eV^] < 3.7xl0-^ 5.4x10"^ < Am| [eV^] < 9.5xl0-^ (2.1) 

at the 3cr confidence level, with Aml^^ = m^^ — m^^ and Am| = m^^ ~ ''^li- 

A minimal, symmetric realization of a model with heavy Majorana neutrinos can be 
obtained by adding to the SM field content one right-handed (singlet) neutrino per family 
UiR, with i — 1,2, 3. The leptonic sector of this minimal model consists of the fields: 




where the obvious labelling, I = (1,2,3) = (e./j.r), will be employed. At temperatures T 
larger than the critical temperature associated with the electroweak phase transition, the 
T-dependent vacuum expectation value (VEV) v{T) of the SM Higgs doublet $ vanishes, 
i.e. ($(T)) = v{T)/\/2 = 0. This is the epoch where a possible leptonic asymmetry created 
by out-of-equilibrium heavy Majorana- neutrino decays can be actively reprocessed into the 
BAU through the equilibrated {B + L)-violating sphaleron interactions. 

At this epoch relevant to leptogenesis, the dynamics of the early Universe is usually 
described by a Lagrangian in the unbroken gauge-symmetric phase of the theory. In this 
unbroken phase, the Lagrangian of the leptonic sector of the model under study may 
conveniently be expressed as 



with 



Aept — -^kin + + -^M , (2.3) 



>Ckin = ^ im^Li + UiRi^UiR + liRi^liRj , (2.4) 
i=i ^ ^ 



i,i=i ^ ' 



(2.5) 



3 



= E [{^iuf {Ms)i^ViK + v,n{Ms)%{viKf\ . (2.6) 
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In the above, £kin, C,y and Cm describe the kinetic terms, the Yukawa sector and the 
Majorana masses of the model, respectively. In addition, $ = 2x2$* is the isospin conjugate 
of the Higgs doublet $, where T2 is the usual Pauli matrix, and the superscript C denotes 
charge conjugation. 

In the unbroken phase of the theory, only the singlet neutrinos are massive. Their 
physical masses can be found by diagonalizing the 3-by-3 singlet Majorana mass matrix 
Ms in ()2.6|) . The matrix Ms is symmetric and in general complex, and can be diagonalized 
by means of a unitary transformation 

U'^MsU = Ms = diag {niN^, mN2, rriNa) , (2-7) 

where f/ is a 3 x 3-dimensional unitary matrix and ?7^7Vi 2 3 denote the 3 physical masses of 
the heavy Majorana neutrinos A''i_2,3, ordered as mAr^ < < "^Afa- Correspondingly, the 
flavour states and (z/j/j)*" are related to the mass eigenstates Ni through 

3 3 
i^iR = Pr E U^JN, , {u^nf = Pl J2 U*,^, , (2.8) 
i=i i=i 

where Pr = (1 + 75)/2 and Pl = (1 — 75)/2. Note that and (i^iR)'" do not transform 
independently of one another under a unitary rotation. In the physical basis, the Yukawa 
leptonic sector reads 

Cy = - (h'^rL^^P^N^ + h[,L,<^Pnh + H.c), (2.9) 

where a four- component chiral representation for all fermionic fields should be understood. 
In fj2.9|) . /i-j is a diagonal positive matrix and h\j is related to h!{f through a bi-unitary 
transformation: If = h'^'^ U, where Vl is a 3-by-3 unitary matrix that transforms the 
left-handed charged leptons to their corresponding mass eigenstates. Our computations of 
the leptonic asymmetries and collision terms relevant to leptogenesis will be based on the 
Lagrangian ()2.9|) . 

Having set the stage, it is now instructive to discuss the possible flavour structure of 
low singlet-scale models with nearly degenerate heavy Majorana neutrinos. Such a class 
of models may be constructed by assuming that lepton-number violation (and possibly 
baryon-number violation) occurs at very high energies at the GUT scale Mqut ~ 10^^- 
10^^ GeV, or even higher close to the Planck scale Mpianck ~ 10^^ GeV through gravitational 
interactions. On the other hand, operators that conserve lepton number are allowed to be 
at the TeV scale. 
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Since our interest is to resonant leptogenesis, the following sufficient and necessary 
conditions under which leptonic asymmetries of order unity can take place have to be 
satisfied by the model under discussion [14]: 



- rriN, 



N,. 



(2.10) 



for a pair of heavy Majorana neutrinos Nij. In ()2.1U|) . T^i are the Ni decay widths, which 
at the tree level are given by 



in 



(2.11) 



In the following, we present a rather generic scenario that minimally realizes the above 
requirements and still has sufficient freedom to describe the neutrino data. Our generic 
scenario is based on the FN mechanism [53]. Specifically, we introduce two FN fields, S 
and E, with opposite U(1)fn charges, i.e. (5fn(S) = — Qfn(S) = +1. Under U(1)fn, the 
following charges for the right-handed neutrinos are assigned: 



Qfn {^ir) 



-1 



(2.12) 



In addition, all other fields, including charged leptons, are singlets under U(1)fn- Then, 
the singlet mass matrix Ms assumes the generic form: 



Ms 



M 



1 

V e 



e 

£ 



\ 



(2.13) 



Mx/M I 



where e = (S)/Mgut and e = (S)/Mgut- In (I2.13j) . M sets up the scale of the leptonic 
symmetry — L^,^ while Mx represents the scale of Lj. violation. It is conceivable that 
these two scales may be different from one another. For the case of our interest, it is 
M ~ 1 TeV, while Mx is considered to be many orders of magnitude larger close to Mqut- 

The FN mechanism also determines the strength of the Yukawa couplings. After 
spontaneous symmetry breaking (SSB), the resulting Dirac-neutrino mass matrix niD has 
the generic form 

/e e 1\ 



V V 



eel 
\e e 1/ 



(2.14) 



^Similar textures of Ms may result from Eg theories [41], where the lepton numbers are approximately 
broken [54]. 
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where is a 3 x 3 matrix containing the neutrino Yukawa couphngs, expressed in the 
positive and diagonal basis of the respective charged-lepton Yukawa couphngs. 



If one assumes that (S) ~ (S) ~ \J M Mqut and Mx ~ Mqut, a rather simple 
pattern for the mass matrices vtiy) and Ms emerges. In this case, the mass spectrum of 
the generic scenario under investigation contains one super-heavy Majorana neutrino, with 
a mass m^vg ~ Mx ~ Mqut, and two nearly degenerate heavy Majorana neutrinos A''i^25 
with rriTVia ~ ^ and a mass difference mx^ — rnx2 ~ ~ M^/Mqut- Since it is 

r^l^ ~ £^^M ~ M^/Mgut, it can be readily seen that one of the crucial conditions for 
resonant leptogenesis in (j2.1U|) . i.e. nix^ — ^atj ~ \ '^n^2i can naturally be satisfied within 
our generic framework. 

In the above exercise, one should bear in mind that the FN mechanism can only 
give rise to an order-of-magnitude estimate of the different entries in the mass matrices 
mo and Ms- Moreover, since our focus will be on the neutrino sector of this minimal 
model of resonant leptogenesis, we will not attempt to explain the complete quark- and 
charged-lepton-mass spectrum of the SM by analyzing all possible solutions through the 
FN mechanism. Such an extensive study is beyond the scope of the present article and 
may be given elsewhere. 

We will now explicitly demonstrate that the mass textures stated in ()2.13j) and ()2.14|1 
can lead to viable light-neutrino scenarios, when the latter are confronted with the present 
solar and atmospheric neutrino data. To further simplify our discussion, we assume that the 
super- heavy neutrino decouples completely from the light-neutrino spectrum. As a result, 
to leading order in the FN parameters e and e, the 3-by-3 light-neutrino mass-matrix 
may be cast into the form: 

/ 2/lii/;,i2 ^11^22 + ^12/^21 /ill/i32 + /i31^12 \ 

hiih22 + hi2h2i 2/2,21/^22 h2ih^2 + ^'-31^22 • (2-15) 

\ hiih^2 + /i3l/il2 /i2l/i32 + /i3l/i22 2/l3i/l32 / 

Here, hij are the neutrino Yukawa couplings in the weak basis described after ()2.14p . Note 
that effects due to the mass degeneracy of the heavy Majorana neutrinos contribute terms 
0{e^e., ee'^) to m'^. As long as e , e < 10"'^, these sub-leading terms do not affect the light- 
neutrino mass spectrum and hence they can be safely neglected. The scenarios which we 
will address numerically in Section 4.3 are compatible with these limits on e and e. 

Let us now present a concrete example by considering the following set of Yukawa 
couplings (given in units oi e = e): 
1 2 

hii = hi2 = -] /i2i = 2; /i22 = 1; /i3i = 1; /i32 = 2 . (2.16) 



m 



2M 
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For our illustrations, we also neglect the existence of possible CP-odd phases in the Yukawa 
couplings. Then, the light-neutrino mass matrix exhibits the structure 



V ^ ^^^^ 



2M 



-4/9 1 0\ 
1 4 5 
V 5 4/ 



(2.17) 



It is not difficult to see that the above light-neutrino mass matrix m'^ can be diagonalized 
by large z/^-z/^ and Ug-i^fi mixing angles, i.e. \0^^i,^\ ~ tt/A and \0,y^uj ~ tt/6. Instead, the 
z/g-i^T mixing angle is estimated to be small, i.e. \9i,^^^\ < 0.1, as is suggested by the CHOOZ 
experiment [47, 55] . Furthermore, the physical light-neutrino masses derived from m'^ are 
approximately given by 

2 — 2 
V EE TTt 

(m,„ m,3) ^ — - (0.04, 1.5, 9) ~ — ^ (0.04, 1.5, 9). (2.18) 

In deriving the last step of (j2.18|) . we have used the fact that \ee\ ~ M/Mqut and rrit ~ 
v/\/2. Observe that although our approach here has been different, the light-neutrino 
masses in ()2.18|) still obey the known seesaw mass relation [44], and scale independently 
of M . In particular, one can easily check that ()2.18j) is compatible with the observed 
light-neutrino mass differences stated in Even though the present example realizes 

a light-neutrino mass spectrum with normal hierarchy, an inverted hierarchy can easily be 
obtained by appropriately rearranging the Yukawa couplings in fl2.16|) . 

Our numerical estimates in Section |3] will rely on neutrino models that make use of 
the generic structures for the matrices Ms and m^), given in ()2.13|) and ()2.14|) . respectively. 



3 Subtraction of RIS and Leptonic Asymmetries 

In this section we wish to address an important issue related to the proper subtraction of 
the so-called real intermediate states (RIS), e.g. heavy Majorana neutrinos iVj, from the 
L- violating 2 — 2 scattering processes. As we will see in Section HI such a subtraction is 
necessary in order to avoid double- counting in the BEs [51] from the already considered 1 — > 
2 decays and 2^1 inverse decays of the unstable heavy Majorana neutrinos. By studying 
the analytic properties of the pole and the residue structures of a resonant L-violating 
scattering amplitude, we are able to identify the resonant part of a 2 — > 2 amplitude that 
contains RIS contributions only. The so-derived resonant amplitude can then be shown 
to exhibit the very same analytic form with the one obtained with an earlier proposed 
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resummation method [14]. Another important result of our considerations is that we can 
define one-loop resummed effective Yukawa couplings that capture all dominant effects of 
heavy Majorana- neutrino mixing and CP violation. 



3.1 Approach to the Subtraction of RIS 

Let us first consider the simple scattering process L$ — N* mediated by a single 

heavy-neutrino exchange A^. This exercise will help us to demonstrate our approach to 
subtracting the RIS part of an amplitude. The more realistic case of resonant leptogenesis 
with two heavy Majorana neutrinos will be discussed later on. To keep things at an intuitive 
level, we assume throughout this section that all particles involved in this process are scalar, 
e.g. scalar neutrinos or sneutrinos Aj that are predicted in supersymmetric theories [56]. 
Nevertheless, we will discuss the complications that may arise in our considerations from 
the spinorial nature of the lepton and heavy neutrino fields. The s-channel contribution to 
the scattering amplitude T(L$ L'"$"'") due to a single A^-exchange reads: 

r.(L* L-*.) = r,(« N') t«(7v- ^ L-tt, , (3.1, 

where the Breit-Wigner-like propagator has been obtained by summing up an infinite series 
of heavy sneutrino self-energies nj^j^(s). The dispersive part of the self-energy Ren^^(s), 
which has been omitted here, can be suppressed by renormalization at the resonant region 
s ~ m~ (see also our discussion below). Instead, its absorptive part Imn^^(m~) = rrif^Tf^ 
is essential to obtain an analytically well-behaved amplitude at s = m~, where F- is the 
total decay width of the heavy sneutrino A^. 

As we will see in Section 0]1, out-of-equilibrium constraints on the heavy-(s)neutrino 
width Fat (F^) imply Fa? <^ (F^ -C m^). In this kinematic regime, the so-called 
pole-dominance or narrow-width approximation constitutes a very accurate approach to 
subtract the RIS part from the squared matrix element |7^(L$ —>■ According to 

the pole-dominance approximation, we have 



(s — m~)2 + rn?~V'^~ 
^ n' N N 



^ 5{s - m'~) e{^) , (3.2) 



(s-m| + im^F^)(s-m|-im^F^) rn^V^ ^ 

where 6{x) and 6{x) are the usual Dirac and step functions, respectively. Notice that 
in ()3.2|) only one residue related to the physical pole at s = m~ — im^T^ is considered 
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by means of the Cauchy theorem. Substituting ()3.2|) into |7^(L$ i.e. after 

squaring p.ip . we can uniquely isolate the RIS part for this process: 

|rRis(^$-^^$^)P = -^6{s-ml)9{^s)\r^{L<l>^N)\'\TB{N^L^<i>^)\'. 

(3.3) 

This last result is fully consistent with the one presented in [51]. As we will discuss below, 
however, the above approach of identifying the proper RIS part of the squared amplitude 
becomes more involved in the presence of two strongly-mixed unstable particles. 

The above derivation of the RIS component of the squared amplitude was based on 
the assumption that the heavy neutrino is a scalar particle. The spinorial nature of 
introduces further complications. It naively violates the factorized form of ()3.3|1 . and 
\Tyiis{L^ —>■ L'^$^)p can no longer be written as a product of a production and a decay 
squared amplitude, i.e. it is not proportional to \Ta{L^ ^"^^ to \Tb{N 

Instead, we find 



niNi. N 



(3.4) 



X J2 ^sis2 Ss3S4 Tr Ta un{p, Si)un{p, S2)Tb TJ un{p, S3)un{p, s^) T} 

^1,2,3,4 

where U]\i{p,s) is the on-shell 4-component spinor, with p"^ = m^, and the trace is un- 
derstood to act on the spinor space. The RIS squared amplitude can be written in the 
factorized form ()3.3j) . only after we perform a Fierz rearrangement of the spinors and have 
integrated over the phase space of the initial and final states in the calculation of the cor- 
responding reduced cross-section (see also our discussion in Appendix |B)) . Then, after the 
phase-space integrations, the only non-vanishing Lorentz structure that survives is of the 
parity-even form: a j5 + 6, where a and b are mass dependent constants. Based on this 
observation, it can be shown that the final result is fully equivalent to ()3.3|) . and amounts 
to substituting into ()3.4|) : 

It is important to note here that the above spin de-correlated subtraction of RIS carrying 
spin can always be carried out independently of the number of the exchanged particles 
in the s-channel, such as heavy neutrinos, provided the resonant amplitude itself can be 
written as a sum of single-pole resonance terms that have the simple factorized form of ()3.ip . 
We will elucidate this point below, while deriving the RIS squared amplitude due to the 
exchange of two heavy Majorana neutrinos. 

Let us therefore turn our attention to the case of two heavy Majorana neutrinos A^i 
and N2. Analytic expressions for unstable particle- mixing effects with three heavy neutrinos 
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Figure 2: Resummed diagrams contributing to the resonant part of the 2^2 scattering 
amplitude of the process L$ — * L'-^(^^ . Depending on the context, L, Ni^2 i^o,y denote scalar 
or fermion particles (see also text). 



are given in Appendix^ Again, we initially assume that the heavy neutrinos A^i and N2 
are scalar particles, i.e. sneutrinos A?"! 2, but we will discuss in Section 3.2 the complications 
originating from their spinorial nature. As is shown in Fig. the s-dependent part Tg of 
the amplitude T[L^ — > may conveniently be expressed as 



%is) = rf An(s)rf + rf Ai2(s)rf + r^A2i(s)rf + ^^A22(s)^^ 



B 



(3.6) 



In ()3.6|) . (^12) represent the vertices $LA'i^2 ('^'^-^'"A''i^2) that include the wave- 
functions of the initial and final states. Analogously with the single heavy-sneutrino case 
described above, Ajj(s) (with i, j = 1, 2) are the corresponding AjAj-propagators obtained 
by resumming an infinite series of heavy sneutrino self-energy graphs [14]: 



An(s) 



A22(S) 



Ai2(s) 



s — m 



Ni 



nn(s) 



n 



12V 



N2 



m~ + 1122(5) 



-1 



A2i(s) 



ni2(s) 



s - my^ + nii(s) 



(s-m|^+nn(s))(s 



(3.7) 



m^^+ 022(5; 



- n?2(5) 



where Iii2{s) = 1121 (s). We assume that the heavy sneutrino self-energies Ilij{s) have 
already been renormalized in the on-shell (OS) scheme, i.e. they satisfy the properties: 

Relliifs) 



Renij(m| 



0, 



lim 







m- 



(3. 



Ni 



where Re indicates that only the dispersive part of the self-energies must be considered. 
Further details on OS renormalization in scalar theories may be found in [50]. 
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To technically facilitate our discussion, it is convenient to introduce the following 
abbreviations: 



Dijis) 
D{s) 



Sij {s 



m- 



Dn{s)D22is) - Du{s)D2i{s) 



(3.9) 



With the above definitions, the resummed A^jAr,-propagators in ()3.7|) can now be expressed 
in the simplified forms: 



An(s) 



A22(S) 



Dis) ' ' D{s) 

In addition, we introduce the quantity 

ZAs] = 



Ai2(s) = A2i(s) 



Puis) 
D{s) 



as 



Dis) ■ 
(3.10) 



(3.11) 



In the OS scheme, with all contributions from unitarity cuts neglected, we obtain the known 
relation for the residues of the diagonal propagators: Z°^(m~ ) = 1. However, Zi{rn~ ) 
is in general complex, but UV finite at order {h'^Y- 

The two complex pole positions s^^ ^ associated with the heavy sneutrinos A^i^2 are 
determined by the equation D{s^^^ = 0, where D{s) is given in ()3.9|) . Since each re- 
summed propagator Ajj(s) given in ()3.10|) contains two complex poles at s = s^^^, it can 
be expanded about ^ as follows: 



An(s) 



A22(S) 

Ai2(s) 
A2i(s) 



D22 


is) 




, D22is) 




+ 




s) 


s^s 


^ ^ D{s) 


JV2 


Zi 


is) 


+ 


Duis) Duis) 


Dui 


s) 


s — 


'n. 


Dnis) Dis) 


Dni 


s) 


Zi 


is) 


+ 


Duis) Z2is) 


D21 


is) 


s — 


'n. 


Duis) s~s^^ 


D,, 


is) 


Z2 


is) 


+ 


D2iis) Zi{s) 


D12 


is) 


s — 




D22is) S-S^^ 


D22 


is) 



+ 



N2 



+ 



+ 



Z,is) Duis) Duis) Z2is) 



s-s^^ D22is) Dii{s) s-s^^ 
Z2is) D2iis) D2iis) Zi{s) 



s-s^^ Dii{s) 



-D22(s) S - S 



+ 



+ 



(3.12) 
(3.13) 
(3.14) 
(3.15) 
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where the elhpses denote off- resonant terms which are non-singular at s = s^^ ^ . Notice that 
we have retained the s-dependent analytic form for the residues in the above complex pole 
expansion by virtue of the Cauchy theorem. Substituting ()3.12j) - ()3.15j) into the s-channel 
amplitude Tg in ()3.(ij) and neglecting off-resonant terms yields 

fs{s) = V,^{s) Vf{s) + V.^is), (3.16) 



with 



= rf(^) - :^r^(^), = r^(^) _ rf . (3.17) 

^22(5) Dii{s} 

Here, it is important to remark that the expressions Vi_2(s) in ()3.17p become identical at 
s = m~ to the resummed decay amplitudes derived in [57], using an LSZ-type reduction 
formalism. Instead, in the present approach, the corresponding resummed decay amplitude 
can be obtained by studying the analytic structure of the residues of the complete resonant 
scattering amplitude, in which the unstable heavy sneutrinos are described as intermediate 
states in the s-channel. The fact that these two approaches lead to identical results provides 
a firm support for the validity of the method developed in [14]. 

The RIS squared amplitude pertinent to the propagation of Ni and N2 can now be 
identified as 

|TRis(L<l>^L^$t)|2 = (3.18) 
\Zi\^7i6+{s-ml) \Z2\^n 5+(s ~ ml) 



with 5+(s-m|^p = 5(s-m|^p e{^). In (jSHD, ^1 (^2) and V/'^ (V/'^) are evaluated 
at s = m~ is = m~ ). Observe that up to higher-order wave-function renormalization 
effects, the RIS squared amplitude (|3.18|) for two unstable particles is very analogous to 
the corresponding one ()3.3p derived for one single resonance. Although we will not address 
this issue in detail here, we simply note that our subtraction approach of isolating the RIS 
part of a squared amplitude can be extended to more than two unstable particles. The key 
observation to be made here is that such a generalization is possible, since the s-channel 
dependent amplitude in the pole- dominance approximation can always be expressed as a 
sum of products of resummed vertices [cf. ()3.17p ] and Breit-Wigner propagators with single 
complex poles [cf. p.ip and ()3.16|) ]. 

One might worry that the subtraction approach described above may not be applicable 
for the case of our interest with overlapping resonances, i.e. for rri]^^ — rriN^ ~ VN-^^/2. 
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However, we should realize that the particles associated with the complex poles, e.g. sat^ j; 
of a transition amplitude have a completely different thermal history, because of their many 
decoherentional collisions with the other particles in the thermal bath. On the other hand, 
the so-called quantum memory effects are expected to play a relevant role only when the 
decay widths FiVia or the mass difference rriN^ — rriNi much smaller than the Hubble 
parameter H governing the expansion rate of the early Universe at T ^ m^r-^ 2 . In the 
former case, one also finds that A^^i 2 are weakly thermalized [20]. Otherwise, our subtraction 
approach not only takes into account the part of the squared amplitude associated with the 
RIS, but also provides a consistent description of the incoherent properties of the heavy 
neutrinos ^. 

3.2 Resummed Effective Yukawa Couplings and 
Leptonic Asymmetries 

Until now in this section, the heavy Majorana neutrinos were mainly treated as scalar par- 
ticles. However, our approach described above for subtracting the RIS from the squared 
amplitude of the process L$ carries over very analogously to a strongly-mixed 

fermionic system, including the spinorial nature of the heavy Majorana neutrinos Ni 
and N2. 

To see the above point, we first introduce an abbreviated form for the one- loop cor- 
rected inverse A^jA'j-propagator matrix: 

Pijir^) = + S,,,(^), (3.19) 

where denote the self-energy transitions Nj{p) Ni{p), renormalized in the OS 

scheme, and p is the 4-momentum of A^i,2- With the aid of these newly-introduced spinorial 
functions (j3.19j) . the resummed A^jA'j-propagators Sij{^) are given by (suppressing the 
argument ^ everywhere) 

Su = (^11- ^12 ^2^2' ^21) '\ ^22 = (^22- ^21 ^12)'', 

512 = — 'S'li ^12 ^22 = ~ ^11 ^12 S22 , 

S21 = -S22 ^21 = -^2^2' ^21 5n , (3.20) 

^For instance, within the context of thermal leptogenesis, the use of a time-integrated CP-asymmetry 
formula, very analogous to the one applied for a coherently oscillating B'^-system, leads to an erroneous 
incorporation of the decoherentional properties of the thermal bath. 
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with Pi/{^) = ^. These expressions coincide with those presented in [14]. 

In analogy to the scalar case, we also introduce the corresponding Z-factors: 

m) = (^Sr^\^)y\ (3.21) 



where the partial derivative d/d ]^ may act on spinorial expressions that depend on |^ and 

Since the heavy-neutrino self-energies Sjj(^) are renormalized in the OS scheme, their 
dispersive parts satisfy the renormalization conditions: 

ReSi,(^) = 0, — ^ Re^iiif, Ui{p) = 0. (3.22) 

Again, neglecting contributions from unitarity cuts, it can be shown [58] that the condi- 
tions ()3.22|1 assure: ^'^^{'f) Ui{p) = Ui{p) and Ui{p) ^°^(/^) = Ui{p). In general, there are 
deviations from this last equality, which result, however, from scheme-dependent, UV-finite 
terms of order (/i^)^. 

The complex pole positions of the resummed A^jAr,-propagators can be determined 
by solving the equations: 

D,{s) = det[Sr,\^)] = 0, (3.23) 

where the determinant is taken over the spinorial components ^. In fact, if det [^ij(^)] 7^ 
for i 7^ j, it is then sufficient to solve one of the two equations: -Di(s) = or D2{s) = 0, 
to find the two complex poles associated with the unstable heavy neutrinos A^^i and N2. 
In particular, exactly as in the scalar case, each resummed A^iA'j-propagator contains two 
complex poles at y/s = y/sj^^^^ = m^/^ " t 

By means of Cauchy's theorem, each resummed A^jA'j-propagator can now be ex- 
panded about the complex poles i/i = i/i^^ ^ follows: 



S - SNi S - SN2 

(^ + m^J ^2 ^ mi ^2i(^ + mjvJ ^1 ^12 ^22' ^ 

022[P) — + + 

S - SN2 S - SNi 



•^To give an example, we note that det — to) = (s — m^)^, with s = p^. Observe that the solutions of 
det — to) = contain a double positive root at ^/s — v^l b. ~ ™j reflecting the fact that the left- and 
right-handed components of a chiral field have the same physical pole as a consequence of CPT invariance 
of the theory. 
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S2M 



'^22 



'12 



S - Siv, S - SN2 



+ 



+ 



(3.24) 



where the exphcit dependence of and pij on j/ is not shown. As before, the elhpses 
indicate that non-singular terms at s = sni2 have been omitted as well as higher-order 
scheme-dependent terms beyond the OS renormalization. 

After inserting the pole-expanded expressions ()3.24|) into the transition amplitude for 
the process L$ we find very analogously to the scalar case: 



TJs] 



ft 



s - SNi 



fi + f 2 



■A (^+mArJ ^2 



with 



A(B) 



A{B) 



^12 ^22 r 



ir^(B) 



A{B) 



(3.25) 



(3.26) 



n 

and yt2^^ =Ti[2^^'^1q- The incoherent subtraction of the RIS from the squared amplitude 
|T(L$ L'"$^)p can be consistently performed, after the de- correlation effect of the 
heavy-neutrino spins discussed in the previous subsection has been taken into account by 
means of (|3.5|) . Then, in the pole- dominance approximation, we obtain for |7Rig(L$ 
L^(^'^)\'^ a formula analogous to ()3.18|) . where | V'i'^2''^'' P is replaced by 



Tr[r^f\:^+m;v,J K2 

Finally, we should stress again that these results are in complete agreement with those 
derived by the LSZ-type resummation approach in [14]. 

In either resummation approach, i.e. the LSZ-type approach or the one followed here, 
the resummed effective amplitudes for the decays Ni2{p) L^, denoted as T/v^j; 
uniquely determined by 



A{Bh 



^2 



— I 

— i 



h'(2 + ^Vt{^) ^-m^,+^S^r(^) 



abs / 



-I abs 
-'21 

-I abs 
^12 



{^)}uN2ip) , 
(3.27) 



where one-loop contributions from the vertices A^^i 2 have also been included. 

In writing ()3.27p . we have implicitly assumed that all Yukawa couplings and masses are 
renormalized in the OS scheme. Then, up to higher-order scheme-dependent terms, only 
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the absorptive parts Vtii^) ^tji^) of the one- loop vertices and self-energies become 
relevant. These are given by [14] 

= A,^P, + Al^^Pn = E (^^Pl + ^^Pr), (3.28) 

i'=i 



IQn IQiT 



= ^^PL = -tE »Pl f(^], (3.29) 



where Aji = A*^ and f{x) = y/x[l — (1 + x) ln(l + l/x)] is the Fukugita-Yanagida one-loop 
function [4]. 

Substituting ()H.28|1 and ()H.29j] into ()3.27|1 and neglecting terms which are formally of 
order {h")'^ and higher yields 

TiN,-^L^) = Ch+)HUiPRUN,, (3.30) 

with 

where the parameters Aij and are defined in ()3.28j) and (j3.29|l . The CP-conjugate 
decay amplitudes T(iVj L*-^$^) can easily be recovered from ()3.30|) by just taking the 
complex-conjugate Yukawa couplings in 1)3.311) . i.e. 

2/1^2* "^A^i {jriN^ AI2 + Al^) 



{h'^-)ii = ^n* + '^Pn ~ 9 9 TZ^A 9 1 

^ ' a zi ^2^^ _ ^2^^ ^ 2z A22 

Notice that as a consequence of CP violation, it is |(/i+)fj| 7^ \(h''_)ii\. In particular, it is 
important to remark that the effective Yukawa couplings hy^ defined in fl3.3ip and ()3.32|) 
contain the dominant part of the one-loop radiative corrections as well as the enhanced 
heavy- neutrino self-energy effects in the kinematic region (m^rj — mArJ -C {rrijq^ + mArJ. 
As we will see below and in Appendix |B1 the leptonic asymmetries and the radiatively- 
corrected collision terms can be expressed in a compact manner in terms of the resummed 
effective Yukawa couplings /i^. 

Let us now turn our attention to the discussion of the leptonic asymmetries ^at^ j. 
The CP-violating quantities 5^^ ^ are defined and easily calculated as 



T{Ni ^ L$) + r(Ar. ^ L^$t) {h\%\)u + {h!'Jh-^^ 



(3.33) 
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L 



Figure 3: e- and e'-types of CP violation in the decays of heavy Majorana neutrinos. 

In (j3.33|) . a matrix notation for the resummed effective Yukawa couplings (h'^)ii should be 
understood, with / = 1, 2, 3 and i = 1,2. As is illustrated in Fig. El it is now interesting 
to discuss the two different types of CP violation contributing to SN^, the so-called e'- and 
e- types of CP violation. If we neglect all self-energy contributions to and 0[{h'^)^] 
CP-conserving terms, we then find the known result for the e'-type CP violation [4]: 



Im (/i^t h 



u\2 



ml, 



with i ^ j. Instead, if all one- loop vertex corrections have been neglected, we obtain a 
simple formula for the e-type CP violation [14,32]: 

^ ih^^ h^h ih^^ h^n {m%^ - m%y + n^vT ' 



where i,j = 1,2 {i ^ j) and the tree-level decay widths Pj^!* of the heavy Majorana 
neutrinos Ni are given in 1)2.111) . In finite-order perturbation theory, the absorptive term 



m 



j^T^^^,'^ that occurs in the last denominator on the RHS of ()3.35p is absent, thereby 
leading to a singular behaviour for e^. in the mass-degenerate limit m^- mi^^.. However, 
the appearance of this regulating absorptive term due to the finite width of the heavy 
Majorana neutrinos should be expected on physical grounds and emerges naturally within 
our resummation approach. Finally, it should be noted that ()3.35|) is valid for a mixing 
system with two heavy Majorana neutrinos only. The generalization of ()3.35p for the case 
of a three-heavy-Majorana- neutrino mixing is more involved and hence has been relegated 
to Appendix |X| 

A non-zero leptonic asymmetry can be created, if and only if the following CP-odd 
quantity does not vanish [14]: 



^cp = Ini Tr 



{K"^)^ h"^ Ml Ms Ml Qf^Y {K^^Y Ms 



^ 0, (3.36) 
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where h'^^ and Ms are the Yukawa-couphng and the singlet neutrino-mass matrices, de- 
fined in the weak basis through ()2.5|) and ()2.6|) . respectively. An important property of the 
CP-odd quantity Acp is that it is invariant under the so-called generalized CP transfor- 
mations [59,60]. In the physical basis, Acp is found to be 

Acp = mNj_mN2{m%^ - m%Jlm {h"^ h^ f^^ , (3.37) 



Acp = ruN-, itin^ {m%^ - m%J Im 



{h''^ K')\2 + 171^2 {m% - m^J (h"^ 



)23 



+ niNi - m%J Im {h"'' h^f^^ , (3.38) 

for a two- and a three-generation heavy-neutrino model, respectively. In general, the total 
number Mcp of all non-trivial CP- violating phases in a model with weak isodoublets 
and neutral singlets is Mcp = nLin^ — 1) [61]. However, after summing over all lepton 
flavours that occur in the flnal states of heavy Major ana-neutrino decays, only one CP- 
violating phase becomes relevant for leptogenesis which can be equivalently represented by 
the rephasing-invariant quantity Acp in ()3.36|) . 

Technically, we may understand this last point as follows. As we will explicitly see in 
Section 4.2, the net source 5ni for generating a non-zero value for the number density rii 
of the lepton number L has the analytical structure: 

driL = J2 PiVi gimN^ , (3.39) 

1=1,2,3 

where PAfi a.s the radiatively-corrected total decay widths of Ni and g{m^^) is an an- 
alytic function of m^Tj that contains Boltzmann-type factors, and whose precise form is 
not important for the present discussion. Since is a physical CP-violating quantity in 
leptogenesis, it should be proportional to Acp. Indeed, it can be straightforwardly checked 
that in ordinary flnite-order perturbation theory, both e- and e'- types of CP violation are 
proportional to Acp, only after the sum over Ni is considered ^. Most remarkably, it has 
been shown for a model with two heavy Majorana neutrinos [14] that even the resummed 
expressions of ^a?. given in ()3.33j] lead to a drii oc A^p; the three-heavy-Majorana-neutrino 
case is discussed in Appendix El This property shows the consistency of our resummation 
approach with respect to generalized CP transformations. 

"^This statement is only valid for a three-generation heavy-neutrino model. For a model with two heavy 
Majorana neutrinos, the e-type CP-violating terms in (SjV; are individually proportional to A^p. 



21 



4 Boltzmann Equations for Resonant Leptogenesis 



Before we solve numerically the relevant BEs, it is useful to give a qualitative discussion of 
the out-of-equilibrium constraints on the parameters of the theory. Specifically, within the 
context of resonant leptogenesis, we demonstrate how moderate departures from the out-of- 
equilibrium condition on the decay rates of the heavy Majorana neutrinos are sufficient to 
significantly lower the singlet scale to the TeV range, without being in confiict with neutrino 
data. Subsequently, we set up the relevant network of BEs, where important contributions 
to scatterings from enhanced heavy-neutrino self-energy graphs and from gauge-mediated 
interactions are included. Solving numerically the BEs for specific scenarios compatible 
with neutrino data, we show that the leptogenesis scale can be lowered up to the TeV scale 
or even lower close to the critical temperature T^, namely up to the scale at which the 
{B + L)-violating sphaleron interactions are still in thermal equilibrium. 



4.1 Out-of-Equilibrium Constraints 

The out-of-thermal equilibrium condition on the heavy Majorana neutrino decays places 
severe limits on the Yukawa couplings of neutrino models [2]. To obtain a qualitative 
understanding of those limits, let us first introduce the parameters 

r(o) p 
K, = (4 1) 

3vel and re 

and H{T) in ()4.ip is the Hubble parameter 



where and Fat. are the tree-level and resummed total decay widths of Ni, respectively. 



H{T) = IMgl/'-^. (4.2) 

J W Planck 

In (jOI), Mpianck = 1-2 X lO^*^ TeV is the Planck mass and ^f* ^ 107 is the number of 
relativistic degrees of freedom of the SM. Obviously, the parameters Ki quantify the de- 
viation of the heavy-neutrino decay rates Fj^:* from the expansion rate of the Universe. 
These parameters should not be much larger than a certain maximum value K^'^^, such 
that leptogenesis can be successfully realized. Even though most analyses conservatively 
assume Kf^^^ ~ 1, much larger values of Kf^'^^ even larger than 1000 can still be tolerated. 
In Section |3J3, we provide a firm support of this observation, after numerically solving the 
corresponding BEs. 
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Let us now consider the out-of-equilibrium constraint Ki < Kf^^^. This is easily 
translated into the upper bound 

{h^^h^)u < 3.5 X 10-1^ {r§v) ■ ^^-^^ 

The very same upper bound can be expressed in terms of new parameters mi, called effective 
neutrino masses in [20], i.e. 

= (^'^^^'^)^^ < lO-3x--eV. (4.4) 
2mNi 

At temperatures T above the electroweak phase transition, i.e. for T > Tc ^ 200 GeV, 
{B + L) -violating interactions mediated by sphalerons are in thermal equilibrium. This in- 
equilibrium condition gives rise to the following relations among the number-to-entropy 
ratios of densities [8,9]: 

YsiT > T,) = ^ Fs-l(T > = Yl{T > T,) , (4.5) 
where we have defined 

Y,^^, (4.6) 
s 

with X — B, L, {B — L), and s is the entropy density. Thus, approximately one-half of the 
lepton-to-entropy density ratio Yl — ul/s gets converted into a baryon-to-entropy density 
ratio Yb — ub/s through the equilibrated {B + L)-violating sphaleron interactions. 

In order to relate the Yb generated at a temperature T = T^, > Tc to t/b measured 
much later at the recombination epoch Tq, we may conveniently assume that there is no 
source or mechanism for significant entropy release while the Universe is cooling down 
from to Tq. Under this plausible assumption of entropy conservation, one can then 
estabhsh the relation Yb{T^.) = Yb{To). Employing now the fact that s(T) = gs{T)n^{T) 
is the entropy density of a plasma with a number Qs of relativistic degrees of freedom at 
temperature T, we arrive at the relation: 

,B = ^ ^ . (4.7) 

For our numerical analysis, we use [2, 20]: gs{To) = 3.91 and gs{T*) — 107, i.e. 
9s{To)/gs{T*) ~ 1/27. If Ki < K^^^, an order-of-magnitude estimate for tjb may be 
obtained by 
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In the above formula, the lepton-to-baryon conversion factor through sphalerons, given 
in ()4.5|) . has also been implemented. It is not difficult to see from ()4.8|) that t]b can 
be around the observed value of 6 x 10~^° if \6Ni\/Ki is of order 10~^. Evidently, CP 
asymmetries of order unity allow for very large values of Ki. In fact, large values of 
Ki ^ 1 lead to a thermally dense plasma, so the required conditions of kinetic equilibrium 
and decoherence of the heavy Majorana neutrinos in the BEs are comfortably satisfied. 
Furthermore, since most of the heavy neutrinos have already decayed at T ^ m^v, ~ 
m^, there are no dilution factors through entropy production at lower temperatures T -C 
rriN [11]. 

It is now straightforward to examine whether the out-of-equilibrium constraints on 
the Yukawa couplings ()4.H|1 . together with the estimate ()4.8|1 for a successful generation 
of the BAU, still allow a light-neutrino sector that can adequately describe the solar and 
atmospheric neutrino data. In the framework of the generic models discussed in Section |21 
their hierarchical light-neutrino mass spectrum requires that Tr (m'^) ^ 0.05 eV, where m'^ 
is given by ()2.17jl . This implies that 



Comparing ()4.3j) and ()4.9j) and assuming no accidental cancellations in the different sums 
of Yukawa couplings, we find that values of Ki larger than ~ 30 are sufficient to successfully 
describe the solar and atmospheric neutrino data, provided Kf^^^ > 100. Most remarkably, 
this result is almost independent of the leptogenesis scale rriN, which can be as low as 
1 TeV. A rigorous demonstration of this observation will be given in Section ^3. 

In our discussion above, we have assumed that finite temperature effects will not 
affect significantly the main results of our analysis. In particular, one may have to worry 
whether the condition for resonant CP violation stated in ()2.10|) will drastically modify 
under the infiuence of thermal effects. Indeed, finite temperature effects on the T = 
masses of the SM particles are significant. Gauge and top-quark Yukawa interactions give 
rise to appreciable thermal masses for the leptons and the Higgs fields [62], i.e. 



where g and g' are the SU(2)i and U(l)y gauge couplings at the running scale T, and 
d = (8M^ + M| + 2m^ + M^)/(8u^). At temperatures T < ttin where leptogenesis becomes 
more operative, thermal mass effects on leptons are roughly one order of magnitude smaller 




(4.9) 



^ (3/ + g") 



J^2 




(4.10) 
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than the heavy Majorana neutrino mass. However, thermal contributions to the mass of 
the Higgs field are more important, and highly depend [63] on the actual value of the T = 
Higgs-boson mass Mh- If the range of Higgs-mass values 115 GeV ^ Mh ~ 190 GeV is 
considered, which is deduced from direct Higgs searches and electroweak precision data, 
one then gets the upper and lower limits: 0.5 ^ M$(T)/T < 0.7. As a result, the effective 
decay widths of the heavy neutrinos T]^.{T = mj^.) will reduce at most by a factor ~ 2 
with respect to Tj^.{T = 0) due to phase-space corrections. 

Instead, thermal effects on heavy Majorana-neutrino masses are very suppressed, as 
they are proportional to the heavy- neutrino Yukawa couplings h^j. Adapting the results 
of [62] to our model, the size of these thermal effects may be computed by 

where Ms{0) is the physical diagonal heavy-neutrino mass matrix defined in fl2.7|l in the 
symmetric phase of the theory at T = 0. At finite temperatures, M|(T) is in general not 
diagonal and therefore needs a T-dependent re-diagonalization. However, from ()4.11|) it 
can be estimated that for nearly degenerate heavy Majorana neutrinos with rriNi ~ ''^n, 
the thermally induced mass splitting is 

mjv,(T) - m^^(T) < 1 Re [(/^^t/,-),,] — . (4.12) 

16 rriN 

This thermally induced mass splitting is comparable to the decay widths Tf^f. of the heavy 
Majorana neutrinos at temperatures T < m^r, at which a net L and B number can in 
principle be created. As a consequence, the conditions for resonant CP violation in ()2.10|) 
are not spoiled by thermal effects in a relevant way. 



4.2 Boltzmann Equations 

Our derivation of the BEs relies on a number of approximations and valid simplifications. 
In particular, we neglect thermal effects on all collision terms which become less significant 
for temperatures T < rriM^ relevant to leptogenesis. In addition, we adopt the Maxwell- 
Boltzmann statistics, which is expected to introduce errors no larger than 20%. Never- 
theless, several crucial improvements, which were neglected in the existing literature, have 
now been implemented in the BEs. To be specific, we include the gauge-mediated collision 
terms that describe processes such as A^jLj^$'''V^ and their crossing-symmetric reactions, 
where collectively denotes the SU(2)l and U(l)y gauge bosons (with a = 1,2,3) 
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and in the unbroken symmetric phase of the SM. Since our main interest is resonant 
leptogenesis, we include heavy-neutrino self-energy enhanced contributions to scatterings 
according to the resummation approach analyzed in Section 121 

Before writing down the relevant set of the BEs, it is useful to establish notation and 
define a number of auxiliary quantities. For this purpose, let us start by reviewing a few 
basic concepts. The number density Ua of particle species a, with ga internal degrees of 
freedom, is given by [51] 

(2vr)= 



na{T) = 9a j 7^ exp [ - (^p2 + m2 - iia{T))/T 



2.^ Hi") ' 

where /!„ is the T-dependent chemical potential and Kn{x) is the rath-order modified Bessel 
function [64]. In our minimal leptogenesis model, the ga factors are: gw^ = 3(?b = 6 and 
g<s> = g<^\ = 2, and for the ith family: g^, = 2, gL, = g^c = 4, gQ^ = ggc = 12, and 
9ui = = 6. If the chemical potential /i^ vanishes, i.e. /i^ = 0, UaiT) automatically 
satisfies an in-equilibrium number-density distribution, which is usually denoted as [11,51] 
rf^{T) and takes on the simple forms in certain limits: 

<\T) = (4.14) 
^ , (m, « T) . 
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Here, we should note that although the condition of thermal equilibrium does not imply by 
itself the vanishing of the chemical potential, the in-equilibrium number density n'^ defined 
for /ia = is a useful quantity in writing the BEs later on. 

In analogy to the formalism introduced in [11], let us define the CP-conserving colli- 
sion term for a generic process X Y and its CP-conjugate one X — ^ y as 

7? = liX^Y) +7(X-.F), (4.15) 

with 

7(X^F) = Jd7rxdTTY{27r)U^^\px-PY) e~P'^/^\M{X ^Y)\\ (4.16) 

In the above, |A^(X —>■ Y)\'^ is the squared matrix element which is summed but not 
averaged over the internal degrees of freedom of the initial and final multiparticle states X 
and Y. In addition, we have used the short-hand notation for the phase-space factors: 
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where Sx = Ui^l is a symmetry factor in case X contains a number nid of identical particles. 
An analogous definition holds for rfvry related to the final multiparticle state Y. Since CPT 
is preserved, the CP-conserving collision term 7y obeys the relation 

lY = Ix- (4.18) 

In addition to the CP-conserving collision term 7y , we may analogously define a CP- 
violating collision term as 

S^x ^ ^(X-^F) - 7(X^F) = -57^, (4.19) 

where the last equality in ()4.19|1 follows from CPT invariance. 

Following [11,51], the BEs for the number densities Ua of all particle species a in a 
given model form a set of coupled first-order differential equations. These coupled differ- 
ential equations can generically be written down as 

drir, 



di 



Ha rix' ny 



(4.20) 



where the sum is over all possible reactions in which the particle a can be annihilated or 
created through a reaction of the form aX' ^ Y ot Y ^ aX' . Special treatment [51] is 
required if a particle species a is unstable and hence allowed to occur as a RIS in a resonant 
process like X — > a — > F. We will discuss below our approach to this problem for the case 
of the unstable heavy Majorana neutrinos Ni. 

In order to reduce the large number of the coupled BEs, we assume that all chemical 
potentials of the Higgs field $, the quarks and the gauge fields are significantly smaller 
than the one associated with the lepton number L. This assumption can be justified from 
a thermal equilibrium analysis of the chemical potentials. For a three-generation model, 
such an analysis yields [9,32]: 

4 15 
/iy = 0, /i$ = — /i^, = --/i^, = (4.21) 

where hl = I]f=i fJ'Li, /^q = J2^=i I^Qi and Hu = Yj\=i f^Ui- Note that the remaining chemical 
potentials fid and fie, albeit comparable to fii, enter the BEs only through sub-dominant 
collision terms, e.g. the 6-quark Yukawa-coupling contribution is smaller at least by a factor 
10^^ with respect to the corresponding one due to t-quarks. 

Employing ()4.2U|) and the definitions ()4.16p and ()4.19p for the collision terms, it is 
now straightforward to derive the BEs for njy^ and that govern the time evolution of the 
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number densities of the heavy Majorana neutrinos and the lepton number, respectively. To 
leading order in the small parameter ni^jn^ ^, we obtain 



"I" eq 



(4.22) 



^ + 3ffn, = - (1 - 5) (ii& - ^7^;^ + s-<ii^ + *7?;' 

+ *7.T' - Slit + *7?'v*' 



2nT 



1l<s> + 27LC$t + 47^t$t + 27q„c + 27^^^ + 27^^ 



+ 2-7^'^" + 2-7^'^ + 2-7^'* + ^ ( ^^^^ + -7^'^ 
+ ^7l<I. + ^7$ty + ^7Ly^ + eq I 7q„c + 7<i,ty 



(4.23) 



In ()4.23p . 7^^t denotes the collision term defined in ()4.19p after the RIS contributions due 
to heavy Majorana neutrinos Ni have been subtracted (see also our discussions in Sectional 
and in Appendix EI) • 

In order to terminate the infinite series of collision terms that could be added in the 
BEs ()4.22p and ()4.23j) . we have developed a systematic expansion in powers of coupling 
constants for all 1 ^ 2 and 2 2 processes ^. More precisely, for all these processes that 
involve only one heavy Majorana neutrino, we have included all collision terms depending 
on the coupling constants as (/i±)^, {h'Q'^g'^, (h'^)'^g''^ and (h'^)'^hl, where h'^ ~ h" are 
the one-loop resummed effective Yukawa couplings calculated in Section |21 In fact, we 
neglected 2^2 scatterings (9[(/;.'^)''] with two external heavy Majorana neutrinos, such as 
NiNj ^ LL and its CP-conjugate part. Instead, we included the collision terms of order 
{h'^)'^ for 2 ^ 2 scatterings where all external particles are massless, e.g. LL ^ 



""The quantity n^^ = 2n^ is the in-equilibrium number density for an individual lepton doublet Li with 

/iL, = [cf. dni]. 

®As we will see in Section0]3, 1^2 processes become important at T < tuni, whereas 2^2 scatterings 
dominate at T > tuni (see also [17,20]). In this respect, 1^3 and 2^3 processes should be regarded as 
higher-order corrections to the corresponding 1 2 and 2^2 processes. 
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An important intermediate step in the derivation of ()4.23|) has been the proper imple- 
mentation of the relations of unitarity and CPT invariance that govern the collision terms 
pertaining to reactions with different number of external particles, i.e. between the RIS- 
subtracted 2 ^-i^ 2 scatterings and 1^2 processes or between the RIS-subtracted 3^2 
processes and 2^2 reactions. More explicitly, on account of unitarity and CPT invari- 
ance, the following perturbative relations for the CP-violating parts of the RIS-subtracted 
collision terms can be established 

7'(g«^-LL$) - 7'(Q^«^L^L^$t) = 5^^i^ + Oim'hl] etc., (4.24) 

where the prime defines an operation of RIS subtraction. Notice that the omission of the 
higher-order terms is fully consistent with our truncated expansion outlined above. In 
this context, we also formally neglected as higher-order effects the CP-conserving RIS- 
subtracted collision terms related to 2 — > 3 scatterings. 

Unlike 2^3 scatterings, 3^2 scatterings, e.g. L^u — > L'-^Q, LL^ Qu'--' etc., 
should not be subtracted, as they have not been counted before. Moreover, CPT invariance 
and unitarity give rise to the following constraints for this set of reactions: 

7(L<l>n ^ L^Q) - 7(L^<l>tn^ ^ LQ^) = 0[ih'')^hl] , 

7(LL<I> ^ Qu^) - ^{L^L^^^ Q^u) = OWYK] etc. (4.25) 

As a consequence, the 3^2 scatterings contribute additional CP-conserving wash-out 
2^2 scattering terms to the BE ()4.23|) . through the resonant exchange of heavy Majorana 
neutrinos Ni. These extra wash-out terms can easily be calculated by applying the narrow- 
width approximation to the 3 — > 2 scatterings, e.g. 

7(LL$ QvP) + 7(L^L^$t _ gC^) = ^75^^ + 0\iyfhl\ etc. (4.26) 

Note that these last relations have already been implemented in the BE ()4.23|1 . Finally, 
since CP violation is predominantly mediated by the resonant exchange of heavy Majorana 
neutrinos A^j, the CP- violating collision terms can be further approximated in terms of the 
CP-conserving ones as follows: 

5: Ni X r Niu'^ r Niu'^ r NiL r NiL , /a r,^\ 
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where the CP asymmetries 6n, are given in ()3.33|) and all CP-conserving collision terms are 
presented in Appendix iBl 

To numerically solve the BE's, we introduce a number of new variables. This will 
also enable us to compare our results with the literature. To this end, we make use of the 
relation between the cosmic time t and the temperature T: 

t = ^ , (4.28) 

2H z = l) ' ^ ' 



where 



T 



(4.29) 



The relation (j4.28|) is valid in the radiation-dominated epoch of the Universe relevant to 
baryogenesis. In addition, we introduce the parameters rja which give the number density 
of a particle species a normalized to the number density of photons, i.e. 



n^{z) 



(4.30) 



with 



n^{z) 



2T^ 



2 m 



1 



7]-2 ^3 



(4.31) 



With the above definitions, the BEs ()4.22|1 and ()4.23j) can be written down in a more 
compact form: 



dr] 



dz 



dz 



H{z = 1) 
1 



N,, 



1 _ !Mi 

^ eq 

■pD (i) I p 

~ Yukawa 



1 -r ± ■ 



Yukawa 



+ r 



Sii) 
Gauge 



H{z = 1) 



1=1 



, f5(i) 

' Gauge 



i=l 



J2 ( r^^'^ + r 



W{i) yWii) 
Yukawa ' Gauge 



S{i) yS{i) 
Yukawa ' Gauge 



+ F: 



AL=2 
Yukawa 



(4.32) 



(4.33) 



where 
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J- Gauge - ^ yiL<l> + 7$ty^ + iLV^ J , 
^ Gauge - 1^7l-I. + ^eq 7$tV^ + IlV^ J , 

J- Yukawa ~ ^ 'Q«C: + IlqC + 7l« + ^Qu^ J ' 

^ Gauge ~ ^ 1^^^* + ^-i-ty^ + 7Ly^ + 2r/^ / ' 

r^u^^i = ;^(7fc%+27|f^t). (4.34) 

In writing the BEs fl4.32|) and ()4.33|) . we used the approximate relations ()4.27p . Hence, all 
CP- violating contributions do only depend on the parameters Sn^^^^. 

Finally, it is worth stressing that the BEs (j4.32p and (|4.33p can be applied without any 
additional restriction to a more general context of thermal leptogenesis, including scenarios 
of non-resonant leptogenesis. Most importantly, we observe that CP-violating scatterings 
provide an additional non-negligible source of CP violation and may lead to an increase in 
the predicted values for rji at T > ttini, as opposed to previous studies where those effects 
were neglected. In the next section, we will present numerical estimates for the generated 
BAU in a few representative light-neutrino models compatible with solar and atmospheric 
neutrino data. 



4.3 Numerical Examples 

In our numerical analysis, we will consider scenarios with two nearly degenerate heavy 
Majorana neutrinos A^^i 2 with masses at the TeV range. The mass of the third heavy 
Majorana neutrino will be taken to be of order 10^^ GeV, so the decoupling of N3 
from the low-energy sector of the theory is natural. For the generic light-neutrino model 
described by ()2.17|) . the neutrino Yukawa couplings {h^)ii^2 are expressed in the mass basis 
as ^ 

ha = {hii £ + /ii2 £ ) , = {hii £ - hi2e^ , (4.35) 

where the Yukawa couplings /i/i^2 are given in ()2.16p . In addition, the FN expansion pa- 
rameters e and e, or equivalently hn and hi2, are assumed to be complex. For the models 
of our interest, \e\ and \e\ are taken to be smaller than 10~^. According to our discussion 
in Section |21 the scenario given by ()4.35p can accommodate the neutrino data, provided 
\e e\ ~ mAT-^/McuT- Therefore, in addition to scenarios with \e\ ~ \e\, we will also present 
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numerical estimates of the BAU for models with \e\ ^ |£:|, where the product \ee\ is fixed 
to uinJMgvt- 

We start our numerical analysis by exhibiting in Fig. |3] numerical values of the lepton 
asymmetry rji and the heavy- neutrino number densities 'r]Ni2 functions of the param- 
eter z = mj\[ /T. Specifically, we have set uin = 1 TeV, xn = -^r^ — 1 = 7.7 x 10"-^°, 
e = 4.3 X 10^^, e = —i4.3 x 10~^. The horizontal dotted line shows the value of t]l needed 
at temperatures T close to the critical temperature Tc = 200 GeV (indicated by the vertical 
dotted line) to produce the observed t]b- In Fig. EJ^a), we have taken r^™ = as the initial 
value of the leptonic asymmetry. As initial conditions for the heavy-neutrino number den- 
sities r^ATi 2 5 "we have considered two possibilities, depending on whether the heavy neutrinos 
are initially in thermal equilibrium, rjf^^^ ~ strongly out-of-equilibrium, ''^jvi 2 ~ ^• 

On the same panel, we also show numerical results obtained if the CP-violating scattering 
terms proportional to 6^^ (ryukawa + ^Giuge) not included in the BE ()4.33|) . The latter 
is the approach followed in the existing literature. We refer to such a numerical compu- 
tation as 'partial'. For the model under discussion, these results are compared with those 
obtained with our 'improved' treatment where those terms are included. We find that the 
two predictions for 77^ at T ^ m^r-^ (2; ^ 1) turn out to come very close together; they 
only differ by ~ 3%. However, at T ~ itln^^ corresponding to 2; ~ 1, the 'improved' and 
'partial' computations may differ even by a factor of 3. Most interestingly, with the new 
CP-violating scattering terms included, the asymptotic value of the leptonic asymmetry rjL 
is attained at a somewhat higher temperature than in the 'partial' approach. 

An important consequence of resonant leptogenesis is that the generated leptonic 
asymmetry rii{T) is independent of the primordial (initial) lepton asymmetry 77™. As can 
be seen from Fig. EJ^b), even if the initial leptonic asymmetry is taken to be maximal 
corresponding to r/™ = 1, this primordial lepton asymmetry gets rapidly erased and the 
predicted rj^is finally set by the resonant leptogenesis mechanism itself for z > 1. Because 
of the independence of r^^ at T = T^, ~ Tc on the initial conditions, in the following we will 
only show numerical values using our 'improved' approach to BEs for initial conditions: 
rfl = and <^ = 1. 

We now present numerical predictions of the BAU for two representative variants of 
our generic model, where we implement the additional constraint, xn = — \ = e^ (with 
e being real), according to our discussion in Section |21 Fig. El shows numerical values of 
'7l, Wi,2 functions oi z = rriN^/T, for two scenarios with uin^ = 1 TeV and rjf = 0, 
rj'^^ ^ = 1: (a) e = e*^e and e = 4.3 x 10"^; (b) e = i^e (with ^ < 1) and \ee\ = 1.85 x lO^^^ 
As before, the adopted values for the FN expansion parameters e and e are chosen to be 
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Figure 4: Numerical estimates of tjl, rj^^^ as functions of z = niM^/T, for a model where 



mjvj = 1 TeV, xm 



-- 1. 



7.7 X 10-1°, e = 4.3 X 10"^ e = -z4.3 x 10"^ and for 



(a) rjj^ = and (b) 77™ = 1. The horizontal dotted line shows the value of rji needed to 
produce the observed t]b- The vertical dotted line corresponds to T = Tc = 200 GeV. "Im- 
proved" and "partial" refer to whether or not the CP-violating scattering terms proportional 

to On. I i v.,V;„„,„ + i Gauge J ^'^^ 



Sn. (r?l„, + r^it ) ^ncluded m the BE CT ). 
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Figure 5: Numerical estimates of rj^, r]Ni2 functions of z = rriN-^/T, for two scenarios 
with niN^ = 1 TeV, xn = — \ = and rjf = 0, rj^^^^ = 1: (a) e = e^^e and 



e = 4.3 X 10 (b) e = i^e and \ee\ = 1.85 x 10 The meaning of the horizontal and 
vertical dotted lines is the same as in Fig. ^ 
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in agreement with neutrino data. In the first variant where \e\ = \e\, we see from|Sfa) that 
CP- violating phases |^| as small as 10"'' are sufficient to account for the observed BAU. In 
this small phase regime, the predicted values for rji and consequently for the BAU scale 
linearly with the CP-violating phase 6. 

The second variant of our generic model realizes a large CP- violating phase {6 = vr/2), 
but introduces an hierarchy between the FN parameters, i.e. e = i^e, with ^ < 1. FromEl^b) 
we observe that a mild hierarchy with values of ^ of order 10~^-10~^ can well explain the 
BAU. As can also be seen from Fig. Efb), the predictions for tjl scale quadratically with 
the hierarchy factor ^ in this scenario. A particularly interesting point that needs be 
emphasized here is that the scenarios with ,^ ^ 2 x 10~^ and 10~^ give rise to out-of- 
equilibrium-departure factors Ki defined in 1)4.11) much larger than 10'^. In particular, we 
find that the simple order-of-magnitude estimate given in ()4.8|) remains valid for values 
at least up to order lO''. 

In Figs.lUandEl we have also displayed the dependence of the heavy-neutrino number 
densities r/ATj 2) temperatures T <C Tc, where the (i? + L)-violating sphaleron interactions 
were turned off. We see that it is t^at-^ j ^ 10~^° at temperatures T ^ 1 GeV. Hence, 
the heavy Majorana neutrinos are under- abundant, much before the epoch of big-bang 
nucleosynthesis [2]. This conclusion holds true for heavy Majorana- neutrino masses as low 
as 0.3 TeV. In this case, the observed BAU can still be generated, with the only difference 
that the reprocessing of leptons into baryons will freeze out at smaller z, i.e. 2; ~ 1. In this 
context, one may wonder whether the leptogenesis scale can be lowered even further, by 
contemplating models where the freeze-out of sphalerons happens at ^ -C 1. Although this 
is in principle possible, one should expect, however, that thermal and non-perturbative 
sphaleron effects will start playing a crucial role for this class of models and hence a 
different approach based on space-time-dependent diffusion equations [65] would be more 
appropriate to reliably address such scenarios. 

We conclude this section by commenting on two points. First, in our numerical 
analysis the gauge-mediated collision terms provide one of the dominant sources of the 
scattering collision terms at T > mTVi- As is illustrated in Fig. IHl the scattering collision 
term Pg^^gg is comparable to Pyikawa which describes the top- Yukawa interactions. Observe 
that Pyiiklwa is negative due to the unphysical RIS-subtracted collision term 7^^t ()4.34|) . 
Nevertheless, it can be shown that the sum Pyiikawa + | ^1=1,2 is in general physically 
meaningful and always positive. Second, we have investigated the theoretical uncertainties 
from the IR mass regulator mjR introduced in Appendix El to deal with IR singularities 
that occur in the reduced cross sections of certain reactions. Specifically, we find that the 
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Figure 6: Functional dependence of the various collision terms contributing to the BEs on 
z = mi\fj^/T. The input parameters are the same as those in Fig. {^a) for 9 = —1.6 x 10^^. 
The corresponding wash-out contributions ^Yu^wa ^'^'^ ^l^^^gg are not shown, as they are 
approximately: Ty^^^^ ^ SFyikLa and T^^i^^ ~ '^^^^nge■ 



uncertainties in the numerical predictions are always less than 10% when miR is varied 
from rriLiT) to m$(T). This last fact provides additional confidence on the stability of the 
numerical results presented here. 



5 Conclusions 



We have studied the scenario of thermal leptogenesis in which the leptonic asymmetries are 
resonantly amplified through the mixing of nearly degenerate heavy Majorana neutrinos 
that have mass differences comparable to their decay widths. We have shown that this 
particularly interesting scenario of baryogenesis, which has been termed here resonant 
leptogenesis, can be realized with heavy Majorana neutrinos even as light as 0.5-1 TeV, 
in complete accordance with the current solar and atmospheric neutrino data. Models 
that might predict nearly degenerate heavy Majorana neutrinos at the TeV and sub- TeV 
scales and lead to light-neutrino mass matrices compatible with neutrino oscillation data 
can be constructed by means of the Froggatt-Nielsen mechanism. Alternatively, specific 
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Eg models [41] in which the lepton number is approximately violated may also naturally 
realize nearly degenerate heavy Majorana neutrinos of TcV mass. 

An important field-theoretic issue we have been addressing in the present article is 
related to the proper subtraction of RIS's from the lepton-number- violating scattering pro- 
cesses. In order to identify the proper RIS contributions, we have examined the analytic 
properties of the pole and residue structures of a resonant L- violating scattering amplitude. 
If the RIS's carry spin as is the case for the heavy Majorana neutrinos, the effect of spin 
dc-corrclation in the squared amplitude should be considered as well. We have shown that 
the present method of extracting the effective decay amplitude from the resonant part of a 
scattering amplitude is fully equivalent to an earlier developed resummation approach [14] 
based on an LSZ-type formalism for the unstable particles. In addition, within the context 
of thermal Icptogcnesis, i.e. with a strongly thermalized bath, all quantum information 
related to the preparation or production mechanism of the initial heavy Majorana- neutrino 
states gets lost. In this case, only the pole positions and residues of a given scattering 
amplitude are invariant under weak-basis transformations. The present approach natu- 
rally embodies these symmetry properties. Hence, it takes consistently into account the 
phenomena of decoherence in the thermal bath of the early Universe, thereby providing a 
natural solution to the so-called initial-state problem. 

Our predictions for the BAU have been obtained after numerically solving the relevant 
network of BE's, where all dominant contributions related to 1 <-> 2 and 2 2 processes 
have been consistently considered. In particular, we have included the enhanced heavy- 
neutrino self-energy effects on scatterings as well as the most important contributions at 
T > rriNi that originate from gauge- mediated collision terms. The self-energy effects on 
scatterings provide new sources of CP violation. As a consequence, the generated BAU at 
T" ~ mjvi could be larger even by a factor of 3, with respect to the one predicted without the 
inclusion of these additional sources of CP violation. Finally, we should stress again that 
our improved BE's are not only valid for a quantitative description of resonant leptogenesis, 
but they can be applied to other thermal leptogenesis scenarios as well, including those with 
hierarchical neutrinos. 

The minimal leptogenesis model under study has the unpleasant feature that it does 
not provide a testable candidate for cold dark matter (CDM). For example, the CDM 
problem could be solved by introducing additional massive sterile neutrinos into the theory, 
but these could not be experimentally observed. Another more attractive solution would 
be to consider supersymmetric versions of the SM with right-handed neutrinos, where the 
lightest supersymmetric particle (LSP), e.g. the lightest neutrahno, is stable because of 
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i?-parity conservation and so it may qualify as CDM. A non-supersymmetric but equally 
appealing alternative would be to consider the invisible Peccei-Quinn (PQ) axion [66-68] as 
CDM [69], thereby solving the known strong CP problem on the same footing. Dedicated 
experiments [70,71] searching for PQ axions can probe this hypothesis. 

Although we have demonstrated that heavy Majorana neutrinos with sub-TeV masses 
can still be responsible for the observed BAU, without being in conflict with neutrino data, 
one may raise the question whether this resonant leptogenesis scenario can give rise to 
further predictions for lepton-fiavour- violating processes [72-74], e.g. for decays /i 67, 
fi — > eee, fi-e conversion in nuclei etc. Here, we should recall that the out-of-equilibrium 
constraints discussed in SectionlHl imply rather suppressed Yukawa- couplings, thus leading 
to unobservably small lepton-flavour-violating phenomena. However, resonant leptogenesis 
allows for significant departures from the out-of-equilibrium constraints, which in turn 
implies significantly less suppressed Yukawa couplings. It would be very interesting to 
study in detail the phenomenological implications of this exciting scenario of resonant 
leptogenesis for low-energy and collider experiments. 
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Note added 

While revising our paper, we became aware of [76] where thermal effects on the collision 
terms, relevant to the domain T > mTVi, were computed. Unlike [76], we have included 
CP-violating contributions to the BE ()4.33|) from 2^2 scatterings which can be several 
orders of magnitude larger than those from 1^2 decays in the above temperature domain, 
at least as is suggested in the T = approximation. In the same context, we have also 
considered the extra 2^2 wash-out terms which originate from the resonant exchange 
of heavy Majorana neutrinos in the 3^2 scatterings. Finally, we reiterate that resonant 
leptogenesis compatible with neutrino data requires sizeable out-of-equilibrium departure 
factors Ki > 30. As a result, 2^2 scatterings thermalize the plasma much faster, so the 
impact of thermal effects on such leptogenesis scenarios for the interesting region T < m^^ 
becomes even less significant with respect to scenarios with ~ 1. 
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A Three-heavy-Majorana-neutrino mixing 



It is interesting to see how the analytic expressions for the resummed A'jA'j-propagators 
Sij{^) given in ()3.20|1 for a model with two heavy Majorana neutrinos generalize to the 
three-heavy-neutrino case. 

Our starting point is the inverse one-loop corrected A^jA'j-propagator matrix: 

srj\^) = s..A^-^N.) + =^.,(^), (A.i) 

where i,j = 1, 2, 3 and is analogously defined to the three- generation case [cf. ()3.19j) ]. 

The inversion of the 3-by-3 matrix- valued matrix in ()A.1|) gives rise to the resummed NiNj- 
propagator matrix Sij{^). In doing so, it proves useful to introduce first the spinorial 
quantities (no sum over repeated indices): 

= mA^) - Pkki^) M) , (A.2) 

where k ^ i and k ^ j, and p^l{^) = [pkkiP)]~^- Making use of ()A.2|1 . the resummed 
A^jA^j-propagators may be expressed as 

= {^'^^-^?^ ^?r^sy' = {^s-^s ^s-'^sy\ (a.4) 

SM) = -s.MP'S P^i-' = -P?^-'p'SsM, (A.5) 

where the dependence of pij on ^ has not been displayed. The remaining entries of Sij{^) 
can easily be obtained by obvious cyclic permutations of the indices i, j, k = 1,2, 3. Exactly 
as in the two-generation case, there are now 3 complex poles for each resummed heavy- 
neutrino propagator. For non-trivial mixing among the heavy neutrinos, the 3 complex 
pole positions can be calculated by ()3.23p . 

The resummed decay amplitudes T(Aj L$) can be calculated by following a line 
of steps similar to those presented in Section 012. For example, we find for T/Vi(A^i —>■ L$) 

Tn, = ui{t,-T, pft' jPf^ - Ta Pft' ) ^nAp) ■ (A.6) 

Here, ri^2,3 contain the proper vertex corrections to ^LNi^2,3- The remaining resummed 
decay amplitudes T/Va^s exhibit an analogous analytic form. Again, only the absorptive part 
of the self-energies and vertices are relevant in the OS renormalization scheme. 
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Exactly as in the two-generation case, we may compute from ()A.6|) the corresponding 
resummed effective Yukawa couphngs h'^ in terms of self-energy and vertex absorptive parts 
Aij and Bu, defined in ()3.28|) and ()3.29|1 and appropriately extended to a three-generation 
model. Neglecting self-energy terms that are formally 0[{h'^Y] in ()A.6j) . we derive the 
resummed effective Yukawa couplings: 



3 



j,k=i 

mN,{mN,Aij + niN^Aji) + Rik \^N,Akj{mN,Aik + mN^Aki) + mNjAjk{mN,Aki + mN^Aik] 



(A.7) 



where 



2 

D,. _ I^^Wi / A . 

mj^^-m'^^+2im%Ajj 



and \eijk\ is the modulus of the usual Levi-Civita anti-symmetric tensor. The respective 
CP-conjugate effective Yukawa couplings (h}'_)ii are easily obtained from ()A.7|) by replacing 
the ordinary Yukawa couplings h'(^ by their complex conjugates. In the decoupling limit of 
m^s "^Afi,2 5 ^A.7jl can be approximated by (jH.Hlj) . Using the resummed effective Yukawa 
couplings (h±)ii derived in this appendix, it is straightforward to compute the leptonic 
asymmetries ^at. in ()3.33|) for the three-generation mixing case. 
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B CP-conserving collision terms 



In this appendix we present analytic expressions for the CP-conserving coUision terms 
required for the numerical solution of the BEs ()4.H2j) and ()4.HHj) . 

For computational convenience and comparison with the literature, the following 
rescaled variables will be used: 

z = —- , X = , ai= { , ar=\ , (B.l) 

where s is the usual Mandelstam variable and miR is an infra-red (IR) mass regulator to 
be discussed below. 

In terms of the resummed effective Yukawa couplings h^^ given in ()3.31|) and ()3.32p . 
the radiatively corrected total decay width Fat. is given by 

= ^ \{K%)u + {h'-^h^M ■ (B.2) 
Using the latter, we may also define the auxiliary parameters Cj as 

c. = f^)' . (B.3) 

We will now employ the formula ()4.16p to calculate the CP-conserving collision terms 
for 1^2 and 2^2 processes that occur in the BEs ()4.32|) and ()4.33|) . These CP- 
conserving collision terms have been defined as 7y in ()4.15|) for a generic process X 
and its CP-conjugate counterpart X . 

For a 1 — > 2 process, e.g. Ni L$ or Ni L'^$^, the corresponding CP-conserving 
collision term 7^4 is found to be 

J (27r)^ EnAp) 



^ A^i V 
TT^ Z 



K,{zy/^), (B.4) 



where qn, = 2 is the number of internal degrees of freedom of Ni, En^{p) = yp'^ + fnj^^, 
and Kn{z) is an nth-order modified Bessel function [64]. 

For 2-^2 processes, the usual definition of the reduced cross-section introduced 
in [11] was used: 

a{s) = 87f<^{s) J d-KY i27cy 6^^\q - py) \M{X ^ Y)f , (B.5) 
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where s = and |A1(X y)p is the squared matrix element summed over all internal 
degrees of freedom of the initial and final states. In ()B.5|) . is the so-called initial phase 
space integral defined as 

^ Jdnx{2nyS^'\px-q). (B.6) 
The above expressions simplify to give the more practically useful equation 

ais) = ^ [dt \M{X^Y)f , (B.7) 

t- 

where t is the usual Mandelstam variable. 

In the calculation of the reduced cross-sections, we face the problem that not all 
of them are IR safe. In processes, such as NiV^ L$, the exchanged particles, e.g. $ 
and L, that occur in the t and u channels are massless. Thus, one finds divergences at 
the t-integration limits t± when performing the phase-space integral in ()B.7|) . A more 
appropriate framework to deal with this problem is finite-temperature field theory, where 
such IR singularities are regulated by the thermal masses of the particles involved in the 
reaction. In fact, although thermal mass effects break the manifest Lorentz invariance, 
they preserve chirality and the gauge symmetries of the theory [75]. In our T = field 
theory calculation of the collision terms, we have regulated the IR divergences by cutting 
off the phase-space integration limits t±, using a universal IR-mass regulator related to 
the thermal masses of the exchanged particles. Evidently, our IR regularization preserves 
chirality and gauge invariance, as is expected from a finite-T calculation. More explicitly, 
for the reduced cross-sections ()B.13|) . ()B.22|) and ()B.23|) which are calculated below, the 
following upper and lower limits of t were used: 

t+ = - "^iR , = - s . (B.8) 

The reduced cross-section ()B.21|) was computed by cutting-off the integral at both t-limits, 
with 

t+ = - "^iR , = m\ + - s . (B.9) 

The IR-mass regulator mjR is chosen to vary between the lepton and Higgs thermal masses, 
rriL^T) and mij,(T), at T ^ m^Vi [cf. fj4.10|) ]. The resulting variations in the predictions of 
the number densities should be considered as theoretical uncertainties due to our T = 
field-theory calculation of the collision terms. 

A simple formula can be found for the CP-conserving collision terms pertinent to 
2 — i> 2 processes, if fjB.5|) and ()B.6|) are inserted into (|4.16|) . In this way, we obtain the 
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Figure 7: AL = 1 interactions between leptons, heavy Majorana neutrinos and quarks. 



expression, 



7# 



m 



(B.IO) 



where Xthr is the kinematic threshold for a given 2 — > 2 process. In addition, the CP- 
conserving reduced cross-section dy in ()B.10|) is defined analogously to 7y in (j4.15p . i.e. 



cr 



Y 

X ! 



(B.ll) 



where the last equality follows from CPT invariance. 

In the following, we present analytic results of CP-conserving reduced cross-sections 
for all 2 — > 2 reactions that contribute to the EEs. As is shown in we start by listing 
the reduced cross-sections involving AL = 1 transitions between leptons, heavy Majorana 
neutrinos and quarks. To leading order in a^, these are given by ^ 



^QuC 



a 



LQC 



3 a,, 



X 



3au {hl%)u + {h''Jh''_) 



In f|Bl2|l and (pl3jl . we have defined 



Tr(/i"t/;,«) 
47r 



^ + ^ln 



2 



Ob Ci^ I (Xf 



(B.12) 



(B.13) 



2M^ 



where and /i" is the up-quark Yukawa-coupling matrix, and rrit is the top-quark mass. 

In addition to the above AL = 1 Higgs-mediated reactions, there are also 2^2 
processes that change the lepton-number by two units, i.e. AL = 2. As is diagrammatically 



^Up to an overall factor, our analytic result in ljB.13|) (specifically the coefficient multiplying the non- 
logarithmic term proportional to —a.i/x) agrees with [11], but differs from the one stated in [20,21]. This 
discrepancy may be traced in the different methods used to regularize the IR singularities. 
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Figure 8: AL = 2 interactions between lepton and Higgs doublets, (a) and (b) correspond 
to the process L$ ^ L*-^<l>"'', (c) anc? (d) correspond to the process LL ^ $^<l>'^. 



presented in Fig |H1 these processes are: L$ ^ and LL and their CP- 

conjugate counterparts. For the first process, particular care is needed to properly subtract 
the RIS's from the reduced cross-section cr^c^^. Xo this end, we first define the Breit-Wigner 
s-channel propagators as 

(B.14) 



X - ai + iy/oiCi 

Then, the modulus square of a RIS subtracted propagator may be determined by 



XI 



TT 



6{x — Oi 







(B.15) 



This subtraction method is in line with the pole- dominance approximation discussed in 
Sectional where the last step of ()B.15|) may be obtained by means of ()3.2|) ®. Note that 
|L)~^(x)p only occurs in the squared amplitude pertaining to an s-channel diagram. With 
the help of the newly-defined quantities in ()B.14jl and ()B.15jl . the properly RIS-subtracted 

^Even though our subtraction approach appears to be similar to the one suggested recently in [76], 
our subtracted RIS propagator squared \D^^{x)\'^ actually differs from [77] by the fact that we use the 
physical spectral representation of H3.2|l for the distribution function S{x — a;), instead of an arbitrary 
regulating function with the same mathematical features. We have checked that the difference of the two 
approaches is numerically insignificant. Therefore, in agreement with earlier remarks in [78] and the recent 
observation made in [76], we also find that the earher RIS-subtraction approaches, see e.g. [11,17,20], tend 
to approximately overestimate the CP-conserving wash-out term in the BE l|4.33|l by a factor of 3/2. 
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reduced cross-section a^c'l.t may be expressed as follows: 
i,j=i I 



(B.16) 



where 



ATx\D'f\ 

X ^/ (^j 



0, 



A 



{St) 



yj dl Cij 

2ttR 



X I Qi rj ^ f X I T 

1 — ^ m ' 

X 



(tt) 



^ Clj Qjj 



271X fa,- — a 



jJ L 



(x + a_,) In ' ^ ~'~ 



(x + Oj) In 



X + Qi 



27ra; 



a,: V a,: 



(B.17) 



The second AL = 2 reaction LL <|)t$t a^ci its CP-conjugate one does not involve 
RIS's and hence can be written down in the shorter form: 



(B.18) 



where 



27r 



Z' aj(x + aj, 



+ 



1 



For i = j, this last expression simplifies to 

1 



In 



(x -|- aj)(x -|- 



(B.19) 



27r 



X 2 Qi ^ / X + tti 

H In I 



X + ai X + 2ai 



(B.20) 



Note that all the expressions for Aij in ()B.17|1 and Bij in ()B.19jl vanish individually in the 
limit X —>■ 0. 

Finally, there are additional AL = 1 reactions that involve the SM gauge bosons 
= B^, and are depicted diagrammatically in Fig. IHl To leading order in ar, the 
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Figure 9: AL = 1 interactions between heavy Majorana neutrinos, gauge bosons, lepton 
doublets and Higgs doublets, (a) and (b) correspond to the process NiVfj, <-> ^L, (c) and 
(d) correspond to the process N^L <-> V^$^ anc? (e) and (f) correspond to the process 



corresponding CP-conserving reduced cross-sections are given by 



^NiVi, _ riygy 
Sttx 



nv9v 



(x + QiY , / X — ai + ttr 

m 



X — ai + 2ar 



^^^^ IGttx^ 



X 



{bx — tti) (oj — x) + 2(x^ -|- xtti — a^) In 



X — ai-\- a., 



(B.21) 



(B.22) 



'LVa 



IQtt x^ 

X{x - tti) 



X — -|- 4ai In 



X — ai-\- 



Or 



(B.23) 



with gv = g',g and ny = 1,3 for = Bj^, W^^, respectively. Notice that up to an 
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overall factor, the analytic expressions for the SU(2)l and U(l)y reduced cross-sections 
are identical. The different overall factors can be evaluated by properly tracing the gauge 
degrees of freedom when a gauge boson couples to an iso-doublet field such as L and $: 
\g^Ti{TaTa) = |<7' for SU(2)i and ^^'^^^(i^l^) = Ig'^ for U(l)y, where ri,2,3 are the 
usual Pauli matrices and I2 is the 2x2 unit matrix. 
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